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The Lieb-Liniger model describes one-dimensional bosons interacting through a repulsive contact potential. 
In this work, we introduce an extended version of this model by replacing the contact potential with a decaying 
exponential. Using the recently developed continuous matrix product states techniques, we explore the ground 
state phase diagram of this model by examining the superfluid and density correlation functions. At weak cou¬ 
pling superfluidity governs the ground state, in a similar way as in the Lieb-Liniger model. However, at strong 
coupling quasi-crystal and super-Tonks-Girardeau regimes are also found, which are not present in the original 
Lieb-Liniger case. Therefore the presence of the exponentially-decaying potential leads to a superfluid/super- 
Tonks-Girardeau/quasi-crystal crossover, when tuning the coupling strength from weak to strong interactions. 
This corresponds to a Luttinger liquid parameter in the range K € (0, oo); in contrast with the Lieb-Liniger 
model, where K € [1, oo), and the screened long-range potential, where K € (0,1]. 


PACS numbers: 03.70.+k, 03.75.Hh, 05.30.-d, 21.60.Fw 

I. INTRODUCTION 

Outstanding developments in the field of cold atoms in op¬ 
tical lattices have opened the path to the experimental design 
and manipulation of many-body quantum states [1]. Specif¬ 
ically, there have been proposals to simulate quantum field 
theories using cold atoms [2, 3]. These proposals are of 
high relevance, as they can access non-perturbative regimes 
of such theories. Regarding the case of long-range interacting 
theories, there have been exciting experimental realizations 
with polar molecules [4] that have enabled the exploration of 
strongly correlated phases not stable with local interaction po¬ 
tentials. 

The theoretical description of such many-body physics is, 
however, challenging. Remarkable techniques have been pro¬ 
posed and applied to a variety of lattice models of strongly 
correlated systems [5-17]. Within this variety of methods, 
tensor networks are a set of ansatze for many-body wave func¬ 
tions proposed to tackle non-perturbative problems in lattice 
models, and have been successfully used to study strongly 
correlated phenomena [18-20]. The continuous counterpart 
of tensor networks has also been subject of recent research. 
More concretely, continuous matrix product states (CMPS) 
have been proposed as a variational ansatz to describe quan¬ 
tum field theories in 1 + 1 dimensions [21]. Another relevant 
proposal to study critical systems is the generalization to the 
continuum [22] of the multiscale entanglement renormaliza¬ 
tion ansatz [23]. 

Applications of CMPS to the study of many-body sys¬ 
tems directly in the continuum include: ground state prop¬ 
erties [21] and excitations [24] of the Lieb-Liniger (LL) 
model, free massive Dirac fermions [25], the A^-flavor Gross- 
Neveu model [25], and two-species of bosonic [26] and 
fermionic [27] systems. In particular for the bosonic case, 
such applications have been focused on continuous mod¬ 
els with contact interactions. Studies of the long-range 
case have been explored before, using other techniques such 
as bosonization [5, 6, 8, 10, 13, 17], numerical computa¬ 
tions [9, 11, 12, 14-16], and perturbation theory [7]. Although 


these approaches have certainly shed light on the physics of 
long-range interacting bosons, most methods have to be re¬ 
stricted to small values of coupling constants or must re¬ 
sort to a discretization procedure. Indeed, interesting phases 
such as superfluidity, Wigner crystal, charge-density wave, 
and Tonks-Girardeau gas have been found and their respec¬ 
tive crossovers/transitions have been discussed. However, to 
the best of our knowledge, there are no studies of bosons with 
long-range interactions, directly in the continuum, using ap¬ 
proaches with no such restrictions on length scales. 

In this work, we use CMPS to study a system of bosons 
with exponentially-decaying interactions. This system can be 
considered as an extension of the exactly solvable LL Hamil¬ 
tonian [28]. Inspired by lattice Hamiltonian models we will 
refer to this bosonic system as the extended Lieb-Liniger 
(ELL) model. We will see below that this extended Hamil¬ 
tonian is a minimal model that captures both the physics of 
the LL Hamiltonian and that of bosons interacting through the 
screened Coulomb potential. When described as a Luttinger 
liquid, we will show that the Luttinger parameter of the ex¬ 
tended Hamiltonian contains values taken by the correspond¬ 
ing parameter in the cases mentioned above. This will imply 
that our extended model displays strongly correlated behav¬ 
ior of both fermionic and bosonic character depending on the 
values of its parameters. 

The outline of this work is the following. In Sec. II we de¬ 
fine the ELL model and show its connection to the LL Hamil¬ 
tonian. Section III is devoted to a scaling analysis of the ELL 
model and its connection to long-range interacting bosons. 
The numerical results based on the variational CMPS will be 
presented and analyzed in Sec. IV. We summarize the main 
results in Sec. V. Technical details of the CMPS are discussed 
in an Appendix. 


II. HAMILTONIAN MODEL 

In this work we will focus on a model of bosonic parti¬ 
cles interacting via an exponentially-decaying density-density 


2 


term in 1 + 1 dimensions. Its Hamiltonian reads 


H = 


1 

2 m 

1/ 


j dx dx'ip^ {x)dxi’{x) — fj, J dxip^ {x)'il^{x) 

dxdy ilj\x)ijj'<{y)ijj{y)'tfj{x), 


( 1 ) 


where tp{x) and 'tjj'’^{x) represent bosonic field operators that 
annihilate and create a particle at point x, respectively. The 
symbol dx stands for the partial derivate with respect to x, 
djdx. The parameter ^ > 0 is the chemical potential, m is 
the mass of the bosons that we set henceforth to 1 / 2 , grj de¬ 
fines the interaction strength of the potential, and p is a char¬ 
acteristic length that controls the range of the interaction. 

Hamiltonian (1) is a particular case of a broader class of 
Hamiltonians that can be written down as 


H = 


1 

2 m 

if 


j dx dxip^ {x)dxi’{x) — fj, J dxip^ {x)'il^{x) 

dxdyw{x — y)'ijj^ {x)ip^ {y)ip{y)’tlj{x), 


( 2 ) 


the fact that in one dimension there is no breaking of contin¬ 
uous symmetries; hence, algebraically-decaying correlations 
are the closest behavior to long-range order. 

Two other instances of Hamiltonian (2), previously ana¬ 
lyzed in the literature, will be relevant to our discussion. The 
first one is the LL model [28], which describes a system of 
bosons interacting through a contact potential weighted by the 
factor g > 0 


wll{x - y) = gS{x - y). (5) 

This model is exactly solvable by Bethe ansatz [28] and can 
be described by Luttinger liquid theory [7, 8 , 30] in the low- 
energy limit. By resorting to dimensional considerations, it 
has been shown that the LL model only depends on the di¬ 
mensionless variable 7 = g/po, which determines the weak 
(7 < 1) and strong (7 ^ 1) coupling limits [28]. The Lut¬ 
tinger parameter as a function of such dimensionless coupling 
7 has been shown to lie in the interval 

if e [1,00), ( 6 ) 


where the interaction potential has been written as an arbitrary 
function w{x — y) of the distance between x and y. The long 
distance behavior of Hamiltonian (2) was studied in Ref. 29 
invoking Luttinger liquid theory [30-32]. A Luttinger liquid 
is characterized by two independent parameters: the velocity 
of the excitations and the Luttinger parameter K. (Luttinger 
liquid theory is equivalent to a conformal field theory of free 
compactified massless scalar bosons [33], with a compactifi- 
cation radius that depends on K.) This parameter controls the 
decay of the correlations and which particular phase governs 
the ground state. 

We define the density-density and superfluid correlations as 


for all values of 7 , resulting in a superfluid ground state. 

The second example is the case of the screened Coulomb 
potential [ 6 , 10 ] 

C 

Wc{z) — ^^2 _|_ J2^q-|-1/2 ’ 

where a > 0 , C defines the strength of the interaction and d 
is a characteristic screening length [ 6 , 10 ] or dimensional cut¬ 
off [5]. Genuine long-range Coulomb interaction is obtained 
by setting a = 0 and d —>■ 0. By performing a bosoniza- 
tion analysis, it has been shown [17] that for a > 0 and as a 
function of C and p. the Luttinger parameter falls in the region 


C{x) 

S{x) 


{p{x)p{Q)) _ 

pI 

(V>t(a;)^( 0 )) 

Po 


if e ( 0 , 1 ]. ( 8 ) 

(3) Moreover, it has been shown [5, 6 ] the existence of a crossover 
from a Luttinger liquid to a Wigner crystal as a —> 0. 
Hamiltonian (1) corresponds to the choice of potential 


where the density is defined as p{x) = tp^{x)tp{x) and 
its expectation value by the ground state is given by po = 
('0^(O)'(/'(O)). In Luttinger liquid theory the asymptotic ex¬ 
pressions for these correlations at long distances, pox ^ 1 , 
are given by [ 8 , 30] 


C{pox) 

S{pox) 


K 1 cos(27rpox) 

27r2 {pqxY ^ {poxY^ 


1 

[poxy/'^^ 


Bo + Bi 


cos(27rpo;r) 

{poxY^ 


(4) 


In these expressions Ai, Bq, and Bi are coefficients that de¬ 
pend on the microscopic parameters in Hamiltonian (2). Here 
we will follow the convention of calling the ground state of 
model (1) superfluid when S{x) decays slower than C{x). 
This will be satisfied if AT > 1/2. Similarly, we will say that 
the ground state has charge order if AT < 1 / 2 , i.e., whenever 
C{x) decays slower than S{x). This convention stems from 


Wexpix - y) = g "'I® (9) 

in Hamiltonian (2). This Hamiltonian is completely speci¬ 
fied by defining the parameter set {g,r],p). In addition to 
the value of the effective interaction strength 7 , potential (9) 
introduces the dimensionless parameter ri/po, which con¬ 
trols its effective range (see Sec. III). In the spirit of lattice 
models of strongly correlated systems, we have dubbed this 
exponentially-decaying interacting system of bosons the ex¬ 
tended Lieb-Liniger model. We will show in this paper that 
by changing the values of (p, p, p), we can obtain a Luttinger 
parameter that covers the interval 

A:e(0,oo). (10) 

Therefore, the ELL model is a minimal model that contains 
the physics describe by both contact potential and screened 
long-range interactions. This is the main finding of our work. 
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FIG. 1. (Color online) Values of the Luttinger parameter K discrim¬ 
inated by the different interactions discussed in the text. For the LL 
model K £ [l,oo). In the case of screened Coulomb interactions 
K G (0,1]. And for the ELL model Tf £ (0,oo). The possible 
phases are quasi-crystal, super-Tonks-Girardeau, and superfluid. The 
Tonks-Girardeau gas corresponds to iC = 1. 


A pictorial representation of this quantity and its relation to 
Eqs. (5) and (7) is shown in Fig. 1. In particular, we will see 
that the LL model results from our extended Hamiltonian on 
the limit p —> oo, and that it also reproduces the long-distance 
physics of the long-range potential for finite 77 . 


III. ANALYTICAL CONSIDERATIONS 

Let us start with a discussion on the physics of the 
exponentially-decaying interactions that will allow us to pre¬ 
dict and interpret the numerical results presented in Sec. IV. 
We begin by noticing that the integrated strength of the po¬ 
tential (9) gives f^^dzwexp(z) = g. In particular, using 
the limit representation 6{z) = lim^_>.oo 2 of the Dirac 

delta function, we see that indeed Hamiltonian (1) reduces to 
the LL model in the rj ^ 00 limit, 

WLhiz) = lim WexpC^)- ( 11 ) 

r]—>-oo 


X^H = [ dx dx'tpHx)dx^{x) — g{\) f dx ipUx)ip(x) + 

2mo J J 


where we have defined the rescaled couplings as 

p(A) = goX, p(A) = poA, fJ.{X) = hqX^, (14) 

and have redefined the ‘bare’ parameters by attaching a sub¬ 
script to them. We have conveniently written down H as X^H 
so we can fix mo without changing the ground state wave 
function. Notice that if we choose X = Pq^ then the di¬ 
mensionless quantity p/po defines the effective range of the 
interaction potential (9). From these scaling equations we 
can deduce that for any given Hq, defined by the parameter 
set {go, rjQ, po), we can obtain a Hamiltonian H\, defined by 
(p(A), 77 (A),/i(A)), by rescaling the couplings as prescribed 
by the relations above. In particular, this implies that if we 
calculate the correlation functions from the ground state of Hq 


It is also enlightening to compare the potential Wexp{z) 
of Hamiltonian (1) to the screened Coulomb case shown in 
Fq. (7). For that, let us Fourier transform Wey^p{z), 

= 7 ^ i + (W’ 

where q stands for the momentum. We can then see that a 
quadratic term Wg^p/g ~ id/v)'^ will be the leading order 

in the long-wavelength limit {q/rj <C 1 ). 

On the other hand, a Fourier analysis shows that the 
screened Coulomb potential can also be written as a quadratic 
function in the long-distance limit provided that a > 1 [ 6 , 10 ]. 
To leading order in q, the interaction potentials of Fqs. (7) 
and (9) are related through g ^ Cld^°‘ and 77 ^ 1/d. Then 
both potentials reproduce the same leading behavior at long 
distance for a > 1. In addition, we notice that exponentially- 
decaying interactions cannot stabilize a Wigner crystal. This 
can be seen by noticing that Fq. (12) cannot produce the lead¬ 
ing order log q, which is a term necessary to stabilize such a 
state [5, 34]. We notice as well that numerical calculations 
have shown that exponential potentials can mimic the general 
behavior of Fq. (7) [35]. 

To explore some further properties of Hamiltonian (1), we 
introduce the scaling transformation: x —>■ x = Xx, for 
the space coordinate x. This change of scale will induce a 
transformation on the field operators, as well, if we want to 
maintain canonical commutation relations: [ij}'''{x),ij}{y)] = 
6{x — y). Explicitly, this transformation is = 

X~^'ijj{Xx), with A = —1/2. We call A the canonical scaling 
dimension of the field Next, we require that Hamilto¬ 

nian (1) remains invariant under the scaling transformation. In 
other words, if H depends on the field 7 / and a set of coupling 
constants labeled by a, then Haii/ix)] = Haltpix)], where a 
is a new set of couplings that in general depend on A. 

Applying the scaling transformation to Fq. (1) we obtain a 
new transformed Hamiltonian, which we write down as 

J dy 7{x)7{y)'4’{y)-i/{x), 

(13) 


we can obtain those for the mapped Hx by rescaling the coor¬ 
dinates of the original correlations [33]. The scaling transfor¬ 
mation also shows that Hamiltonian (1) has actually not three 
but two independent parameters, since they can be related by 
the scaling factor A. For instance, we can eliminate 77 ’s equa¬ 
tion and write down 

g{X) = ^ 77 (A), 77 (A) = ^ g{Xf. (15) 

Vo Vo 

As discussed at the beginning of this Section, the LL Hamil¬ 
tonian is contained as a limiting case of the ELL model. This 
limit corresponds to 77 (A) —>■ 00 which in turn implies for the 
bare couplings that 770 —>■ 00, or 770 —)• 0 and po 0 , such 
that g{X) and p{X) remain finite. In the space of parameters 
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FIG. 2. (Color online) Entanglement entropy as a function of the 
bond dimension D for 7 = 0.83 and T]/po = 3.29. We have fitted 
the data to the scaling relation Sp = ("InD+ 0(1/InD), where 
C = [sJVlIc + 1)“^ and c is the central charge [38]. 



FIG. 3. (Color online) Superfluid correlation function as a function 
of distance for the LL model {t] — >■ 00 ) with 7 = 0.75, and several 
values of D. Inset: Corresponding finite-entanglement scaling of the 
extracted K. 


defined by ( 5 , rj, p) the LL case is constrained to the plane 
p(A) = 00 where g{X) and /i(A) can take arbitrary values. 
As we shall see below, our numerical results indicate that for 
large, albeit finite, values of 77 (A) the ground state of Hamilto¬ 
nian (1) behaves in a similar way to that of the LL model, for 
the same values of ( 5 (A), /i(A)). In the opposite limit, for val¬ 
ues of 77 (A) ~ 1, the physics of the ELL model is not related 
to that of the standard LL Hamiltonian. 

Similar conclusions on the scaling properties of Hamilto¬ 
nian ( 1 ) can be drawn by exploring the scaling of the di¬ 
mensionless couplings 7 = go/po, po/po, and po/po- By 
introducing the transformation under change of scale of the 
density p(A) = po/X, the rescaled couplings in this case 
transform as 7 (A) = 7 oA^, p{X)/p{X) = (? 7 o/po)A^, and 
p{X)/p{X)‘^ = (po/Po)^'^- The resulting set of independent 
equations now read 


7(A) 


70 f t?(A) \ 

Vo/po \pWJ 


m(a) 

p{xy 


Po f pWV 
\pW) 


(16) 


The ELL model now defined by the couplings ( 7 , 77 /p, p/p^) 
will correspond to the LL model when 77 (A) /p{X) 00 . And 
as before, for p{X)/p{X) ~ 1, the ELL Hamiltonian will give 
rise to different phenomena than that of the LL model. 


IV. NUMERICAL RESULTS 

We now discuss the numerical results obtained for the ELL 
model using the CMPS by Verstraete and Cirac [21, 36], 
along with the time-dependent variational principle proposed 
by Haegeman et al. [25, 36, 37] as minimization method. 
The CMPS method produces an approximation to the ground 
state wave function. Prom this approximation we can com¬ 
pute quantities such as the ground state energy and particle 
densities as well as correlation [21, 25, 36] and spectral [24] 
functions. Por our translationallly invariant system, this varia¬ 
tional ansatz is parametrized by two D x D matrices Q and R 


(see Appendix for details). These matrices contain the varia¬ 
tional parameters of the ground state wave function. The bond 
dimension D is a refining parameter that permits to control 
the accuracy of the resulting wave function (including expec¬ 
tation values and correlation functions). The computational 
cost grows as 0{D^). 

We have studied ground states for values of the bond dimen¬ 
sion in the range D = 4 — 32. As a first check of the accuracy 
of the CMPS wave function in describing the ELL model, we 
have calculated the entanglement entropy in order to extract 
the value of the central charge. A typical result for the entan¬ 
glement entropy as a function of D is shown in Pig. 2. Using 
the finite-scaling entanglement formula for the entanglement 
entropy, proposed in Ref. 38, it is possible to extract an esti¬ 
mated value of the central charge of the conformal field theory 
underlying the ELL model. Such theory corresponds to the 
Luttinger liquid which possesses a central charge c = 1. The 
extracted value from the numerics c « 0.95 gives rise to an 
error of 5% and compares fairly well with the expected result. 
We interpret this result as a confirmation that the CMPS gives 
an accurate approximation to the ground state wave function 
of the ELL model. In particular, this statement implies that 
the CMPS is capable of describing critical theories such as 
the Luttinger liquid, which is the low-energy effective theory 
of the ELL model (see Sec. II). 

Within the range of values of D studied, we have observed 
convergence of the quantities calculated in this paper, includ¬ 
ing the Luttinger parameter K. Pigure 3 shows an example 
that correlation functions converge up to distances around 50 
times the interparticle spacing, x ~ 50po- This plot also ex¬ 
hibits an increasingly large power-law region as D increases. 
Likewise, K versus D shows a systematic convergence. K 
has been extracted by fitting our numerical results of the cor¬ 
relation functions to Eq. (4). The fitting variables are K, Ai, 
Bq, and Bi [39]. The extrapolation of K{D) to D ^ 00 
produces an estimate of the error of K for a given D. In the 
example of Pig. 3, assuming that K{D) is quadratic in 1/D, 
this error is 5% for D = 24. Similar errors are obtained in the 
rest of results discussed in this paper. These results have been 
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FIG. 4. (Color online) Weak coupling limit: 7 = 0.75. Superfluid 
S = S{pox) (full) and density C = C{pox) (dashed) correlations as 
a function of the distance for several values of rj/po- The dash-dotted 
lines are results for the LL case, see Eq. (5). The continuous and 
dashed lines correspond to the ELL Hamiltonian (1). A superfluid 
ground state is found in this regime. 


EIG. 5. (Color online) Luttinger parameter K as a function of the 
potential effective range p/po at weak coupling. The values of K 
extracted from the fittings to the data in Eig. 4 for the superfluid 
(charge density) correlations are labelled as SF (CC). For the LL 
model {p —>■ 00 ) the label LL corresponds to the values of K ex¬ 
tracted from bosonization formulae (4) [ 8 ]. 


calculated with D = 24 unless otherwise stated. 

In the following, we will see that exponentially-decaying 
interactions lead to Luttinger parameters in the range of 
Eq. (10), thus containing the cases of screened Coulomb 
[Eq. ( 8 )] and contact [Eq. ( 6 )] interactions. Having K 
spanning such range will lead to crossovers from superfluid 
to super-Tonks-Girardeau to quasi-crystal states. Here we 
will refer to the super-Tonks-Girardeau regime as a state 
of suppressed superfluidity which can be described as spin¬ 
less fermions (the Tonks-Girardeau gas) interacting repul¬ 
sively [7, 14]. 

We will of focus on two relevant regimes: weak coupling, 
where 7 < 1 , and the strong coupling limit, for which 7 1 . 

We have fixed the chemical potential to /i = 0.5 throughout. 
Results for other values of the couplings are connected via a 
scaling transformation, as discussed in Sec. III. Eor each lim¬ 
iting case, we have varied the value of the characteristic length 
of the potential p and follow the evolution of the correlation 
functions. 


A. Weak Coupling 

Let us start by analyzing the weak coupling limit for 7 = 
0.75. Eigure 4 shows the results for the density C{x) and 
superfluid S{x) correlation functions as a function of the dis¬ 
tance scaled to the density pq. As can be seen, the superfluid 
correlations closely resemble those of the the LL case (shown 
as the dash-dotted black line). This suggests that the values of 
the Luttinger parameter of the extended model will be presum¬ 
ably similar to the LL model, and only differences will appear 
at short distances, pox < 1. The density correlations confirm 
these findings. Indeed, for pqx ^ 1 the long distance behav¬ 
ior of the extended and standard LL models is closely related. 
Again differences are found at short distances p^x ~ 1, i.e. at 
energies not reachable with the theory of Luttinger liquid. 


The resulting values of K at weak coupling for the data of 
Eig. 4 are shown in Eig. 5. Our results for the LL case {p -p- 
00 ) compare fairly well to the values of K predicted using 
bosonization, giving an error of ^ 5%. As was discussed 
for Eig. 4, we observe an overall superfluid state in the ELL 
Hamiltonian at weak coupling; meaning K > 1 according to 
Eq. (4). Eigure 5 shows that AT > 1 for the ELL Hamiltonian, 
leading to S{x) decaying the slowest and hence to a superfluid 
ground state. Asp/pQ grows, the value of K is closely related 
to that of the LL model, as expected. 

It is possible to show what the possible values of K are 
for the ELL model in the weak coupling limit, see Eq. (10). 
Starting at 77 — 00 , we know that the exponentially-decaying 
interactions (9) contain the LL model. On the other hand, 
for finite but large p/po we have observed as well superfluid 
behavior {K > 1), and as p /pq decreases the prefactor gp in 
Hamiltonian (1) will be small, leading to quasi-free bosons, 
which imply K —i' 00 . Consequently we can expect that, at 
least in the weak coupling limit, the ELL model will have a 
Luttinger parameter in the region AT G [l,c)o). Notice that this 
interval is the same for AT in the LL model show in Eq. ( 6 ), 
see also Eig. 1. 


B. Strong Coupling 

The strong coupling limit results of the correlations are 
shown in Eig. 6 . Eirstly, for the LL case where p 00 
(dash-dotted line), the dominant correlations are those of a 
superfluid ground state, similarly to weak coupling. The only 
remarkable difference is that for 7^1 the charge correla¬ 
tions display Eriedel oscillations characteristic of the Tonks- 
Girardeau regime [9, 40], where the bosonic system maps to 
free spinless fermions [7, 14]. Secondly, for the extended 
model at large p/po, superfluidity is suppressed although still 
remains as the dominant fluctuation at long distances. How- 
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FIG. 6 . (Color online) Strong coupling limit: 7 = 184. The notation 
is the same as in Fig. 4. In this limit the superfluid correlations are 
strongly suppressed whereas the charge correlation function displays 
an increase indicating a crossover from the super-Tonks-Girardeau to 
a quasi-crystal state. 


ever, charge correlations become increasingly large in the 
short distance. Notice that for large but finite ry/po, super¬ 
fluidity is greatly suppressed compared to the LL model. An 
overall increase of Friedel oscillations is also observed in the 
charge sector. 

Further decreasing of p/po leads to an almost complete 
suppression of superfluid correlations and charge fluctuations 
thus govern at short and long distances. In the range where 
charge correlations decay slower than superfluid ones, we ob¬ 
serve the appearance of a definite wavevector Q that modu¬ 
lates the density fluctuations. As expected from the correla¬ 
tion functions (4), this wavevector is set by po. The average 
density sets a length scale a = 1 /po and we can associate a 
wavevector to it as Q = 27r/a = 27rpo. The appearance of 
this wavevector in C{x) signals the establishment of charge 
order. 

Figure 7 displays the results of the Luttinger parameter at 
strong coupling. As can be seen, the value of K is restricted to 
K < 1, which indicates that the ensuing physics is not related 
to the LL model. As p/po in varied, K increases and presum¬ 
ably when T] ^ 00 the results of the LL model are recovered. 
The estimated error of K for the LL model is around 3%. By 
inspecting Eq. (4) a crossover from a state with suppressed su¬ 
perfluidity to a quasi-crystal state is obtained when K = 1/2. 
Accordingly, we expect that upon decreasing p/po (i-e., away 
from the LL limit) the tail of the exponential interaction will 
dominate, leading to the formation of the ordered state for 
K < 1/2. For K > 1/2 the bosonic system is described by 
strongly interacting repulsive spinless fermions; this is a sup¬ 
pressed superfluid state or the super-Tonks-Girardeau regime. 

The overall picture is then that, for a fixed value of 7 and 
by varying p/po, a crossover from a super-Tonks-Girardeau 
to a charge-ordered state is observed. This ordered state is in 
fact a quasi-crystal, which is the closest state that resembles 
a Wigner crystal within Luttinger liquid theory. The presence 
of the super-Tonks-Girardeau and quasi-crystal states, which 
are not present in the LL Hamiltonian, is directly related to 



FIG. 7. (Color online) Luttinger parameter if as a function of the 
potential effective range p/po at strong coupling. The values of K 
extracted from the fittings to the data in Fig. 6 . The rest of the nota¬ 
tion is the same as in Fig. 5. The crossover value at K — 1/2 from 
the super-Tonks-Girardeau limit to the quasi-crystal state is shown as 
a dashed line. See text for details. 


the exponentially-decaying potential (9). 

We have performed similar calculations for larger values 
of 7 (not shown here). The general tendency is similar to 
that shown in Fig. 6 . However, the range of values of K is 
even smaller than the values reported in Fig. 7. This sug¬ 
gests that in the strong coupling limit of the ELL Hamiltonian, 
the Luttinger parameter can range in the interval K G (0,1]. 
This assumption is supported by the discussion of Sec. III. 
There we have shown that Hamiltonian (1), with specific val¬ 
ues of (p,p), can describe screened Coulomb interactions in 
the long-wavelength limit (see Eig. 1). In addition, bosoniza- 
tion results on systems interacting through potential (7) have 
shown that the compactiflcation radius lies in the same region 
as AT in the ELL model for 7 1. Notice that such values of 

AT at strong coupling match those of Eq. ( 8 ). 


C. Luttinger Parameter 

The analysis and conclusions drawn in Subs. IV A and IV B 
can be further substantiated by calculating the Luttinger pa¬ 
rameter for arbitrary values of {g,r]), with p = 0.5. Re¬ 
sults for other coupling values are connected via the scaling 
transformation discussed in Sec. III. In Eig. 8 we show the 
Luttinger parameter, for the ELL model, versus 7 for several 
color-coded intervals of p/po- The data confirms our expecta¬ 
tions that AT G (0, 00 ), showing that Hamiltonian (1) contains 
as limiting cases the LL model and the screened Coulomb po¬ 
tential. 

As discussed for the weak coupling (7 <C 1) case in 
Subs. IV A, the behavior of AT for the ELL model matches 
that of the weak-coupling regime of the standard LL (full line 
in Eig. 8 ) model. This indicates that the effect of the poten¬ 
tial range does not greatly affect the physics of the extended 
model for 7 ^ 1. In this regime, superfluid correlations gov¬ 
ern the ground state, i.e., K > 1. 
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FIG. 8. (Color online) Luttinger parameter K versus the dimension¬ 
less interaction 7 = g/ po for color-coded intervals of r]/ po, for the 
ELL model. Full (open) symbols correspond to D = 18 (24). The 
values of the effective range of the potential span r]/ po € [0.3, 23]. 
The full line is the bosonization result for the standard LL model, 
obtained with Luttinger liquid theory [ 8 ]. The dashed lines cor¬ 
respond to the superfluid/super-Tonks-Girardeau and super-Tonks- 
Girardeau/quasi-crystal crossovers at K = 1 and K = 1 /2, respec¬ 
tively. 

At intermediate coupling, 7^1, the exponentially- 
decaying potential starts changing the level of correlations, 
and a departure from the LL result is observed, depending on 
the value of ry/po- Indeed for large p/po ^ LL-like trend is 
still seen, in accord with the results discussed in Subs. IV A. 
On the other hand, for smaller p/po the Luttinger parame¬ 
ter crosses the Tonks-Girardeau point K = 1 into the super- 
Tonks-Girardeau state, where superfluidity is suppressed (see 
Subs. IV B). 

Finally, for strong interactions, where 7 1, a crossover 

from the super-Tonks-Girardeau regime to the quasi-crystal 
state is seen as a function of both fixed p/po and increas¬ 
ing 7 and fixed 7 and decreasing p/po. The crossover line 
is defined by K = 1/2 (see Subs. IV B). Figure 8 shows 
that for large p/po the strong coupling limit of the extended 
model tends to closely follow that of the LL model. As p/po 
decreases an overall superfluid/super-Tonks-Girardeau/quasi- 
crystal crossover is observed. 

Based on the previous discussions, we conclude that for 
arbitrary values (g, p, p) the low-energy physics of the ELL 
Hamiltonian is described by Luttinger liquid theory with a 
Luttinger parameter lying in the range (10). This result is in 
high contrast with the cases of contact [Eq. (6)] and of power- 
law [Eq. (8)] interactions (see Fig. 1). This entails that by 
tuning (g, p, p) the decaying exponential can be short-ranged 
enough to describe LL physics or power-lawed enough to 
obtain an analog behavior to that observed for the screened 
Coulomb interaction. 


V. CONCLUSIONS 

Let us summarize the main results presented in this paper. 
Firstly, we have introduced a model for bosons in 1 -F1 dimen¬ 


sions interacting through an exponentially-decaying potential. 
Secondly, by employing well-established scaling transforma¬ 
tions we have shown that this ELL model contains, in some 
limiting cases, both the standard LL model and the long- 
wavelength limit of the screened Coulomb potential. This dis¬ 
cussion allowed us to make some predictions on the phases of 
exponentially-decaying interacting bosons such as the pres¬ 
ence of quasi-crystal and super-Tonks-Girardeau states. Such 
states are not present in the original LL model, for which the 
ground state is always superfluid. Thirdly, making use of the 
recently developed CMPS techniques we have explored the 
ground state phase diagram of our system of exponentially- 
interacting bosons. 

By calculating the superfluid and density correlation func¬ 
tions, we have shown that at weak coupling superfluidity gov¬ 
erns the ground state in much the same way as in the LL 
model. At strong coupling, however, superfluidity is strongly 
suppressed with a simultaneous increase of density correla¬ 
tions, signaling the emergence of the super-Tonks-Girardeau 
state. Upon increasing the interaction density correlations 
dominate the ground state and a quasi-crystal state is sta¬ 
bilized. As noticed above, these additional phases are not 
present in the original LL Hamiltonian. Hence, the decay¬ 
ing exponential potential induces a crossover from superfluid 
to super-Tonks-Girardeau to quasi-crystal states, when the in¬ 
teraction strength is varied from weak to strong coupling. Ei- 
nally, resorting to Luttinger liquid theory, we have shown that 
the value of the Luttinger parameter ranges in the interval 
K G (0, 00 ); thus differing from the values taken by the corre¬ 
sponding quantity of both the LL model, where K G (1, 00 ), 
and the screened Coulomb case, for which K G (0,1). 
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Appendix: Uniform Continuous Matrix Product States for 
Exponential Interactions 

In this paper we numerically investigate translationally in¬ 
variant held theories of interacting bosons in 1 + 1 dimen¬ 
sions, using the technique of continuous matrix product states 
(CMPS) of Verstraete and Cirac [21]. In this appendix we give 
a more detailed discussion of the formalism and its implemen¬ 
tation. 
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Similar to the case of MPS [18, 20] on a lattice, the method 
makes an ansatz for a ground state wave function |4') in terms 
of a set of continuous, matrix-valued functions Q{x) and R{x) 
on an interval x € [—L/2, L/2\. 

Here, Q{x) and R{x) sxe D x D matrices for every point 
x (comprising the variational space), vi and Vr are boundary 
vectors at x = ±L/2 which incorporate boundary conditions, 
Ve is the path-ordered exponential, 1 and (x) are identity 
and creation operators acting at position x in space, and | 0 ) 
is the vacuum defined by '0(x)|O) = 0. D is called the bond 
dimension of the CMPS. When going to the thermodynamic 
limit L —>■ oo (see below), vi and Vr will drop out of any 
equations and can be neglected. For later reference we define 

ry 

U{x,y) =Vexp / dx (5(x) 0 1-f i?(x) 0 '0’^(x). (A.2) 
J X 

The CMPS Eq. (A.l) can be considered to be the limit of a 
certain type of lattice MPS: consider a discretization of the in¬ 
terval [—L/2, L/2] into N equidistant points x„, separated by 
e. One can show [21] that when expanding the path-ordered 
exponential, defining c\ = y/eip’’'{xn) and collecting orders 
of e, the resulting expression is equivalent to the one obtained 
from an MPS 

!</) = ^ • • • (4)*"|0) (A.3) 

{ir.} 


with matrices A*" restricted to the form 


for 1$) is given by locally varying Q{x) and i?(x), and taking 
a superposition of all these variations [37]: 

nL/2 

1$)=/ dxU{-L/2,x) 

J-Lj2 

X (C(x) C) 1 -f W(x) 0 tp^)U{x,L/2)\0). (A. 6 ) 

V(x} and W{x) are the variations of Q(x) and R(x), respec¬ 
tively. When added to |T'), the resulting state is again a CMPS 
of bond dimension D. Such vectors are also referred to as 
tangent vectors [36, 37]. The optimal V*{x) and W*{x) are 
determined by minimization: 

{V\W*} = (A.7) 

For translationally invariant systems in the thermodynamic 
limit (L —oo), which is what we consider in the following, 
Q{x), R{x), W{x), and V{x) can be chosen to be indepen¬ 
dent of X. 

A gauge transformation [18, 20, 36] for a CMPS is a 
transformation on {Q,R) which leaves |T') invariant. It in¬ 
duces a redundancy in the tangent space, because it im¬ 
plies [36, 37, 41] the existence of certain choices of non¬ 
zero Vo and Wo such that the resulting tangent vector |$) 
to the state jT') is zero: |<l>[Vo, Wg]) = 0. Such undesirable 
variations of j'k) can be excluded by choosing a particular 
parametrization of V and W (other choices are possible [36]): 


w = 


(A.8) 


A*" ® — 1 -f e Q(xn) 

Agfc (A.4) 

= ^R\xn). 

Like for lattice MPS, the goal is to approximate the ground 
state wave function of a Hamiltonian H = J dx h{x) in terms 
of a CMPS. Time evolution is done using the time depen¬ 
dent variational principle, proposed by Haegeman et al, for 
CMPS [37]. In this method, time evolution is carried out by 
constructing ^ I'I'(t)) = and using it to update 

the wave function |T'): 

|T'(r-f dr)) = |T'(t)) — driJjT'). (A.5) 

The time dependent variational principle teaches us that the 
best optimal approximation to iT |T') is given by a tangent vec¬ 
tor (see below) belonging to the tangent space to variational 
manifold of the CMPS. This implies a projection of iT|4') 
onto the tangent space, at which point the procedure becomes 
approximate (see Ref. 37 for details). 

The goal is thus to find a tangent vector | $) of given hxed 
bond dimension D which optimally approximates 7714/), such 
that |T') — (7 t| 4>) is again a CMPS of bond dimension D. Usu¬ 
ally 77 is a sum of local operators and thus a reasonable ansatz 


I and r are the left and right reduced density matrices [36], 
obtained from solving the equations 

^ = {l\T = {1\{Q + 1 ® Q + R® R) 

= iq + qH + rHr = o 

d\ ) 

--£7 =T|r) = (QOl-f l(g)Q-f7?(g)7?) |r) 

= Qr + rQ^ + RrR^ = 0. 

T is called the transfer operator, and acts as a superoperator 
on the vectors (Z| and jr). I is in this respect a reordering of 
the vector (^| into a matrix. We use the convention {l\A®B = 
BHA and A® B\r) = ArB^ . To order e, the operator e'^'^ 
equals the MPS transfer matrix E — [19]. 

For our simulations we use the gauge freedom [2 1 ] to hx the 
gauge of [47) such that I = 1. Eq. (A.8) enforces (<i>|T') = 0; 
furthermore, we have 

($1$) = (5(0)tr(ry^). (A.IO) 

We will consider the Hamiltonian (1). Using parametriza¬ 
tion (A.8), ($|77|T') can be evaluated to [24, 41] 
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= 5(0) 


ai 


2m 


[Q,R] 0 [Q,R] - fiRi^ R + R® RC[w]{-T)R'Si R (-T)p^ (l 0 F + i?0 W") |r) 


{1\{R 0 R)£[w]{-T) {1<»V + R<»W) £[w]{-T){R ® R)\r) + R] ® ([0, + [t^, R]) 

- fiR^W + {R(£) W)£[w]{-T){R ®R) + {R(g) R)£[w]{-T){R 0 T^)) |r) 


(A.ll) 


where £[ui](—T) = dzw{z) corresponds to the 

Laplace transform of the interaction potential w{z) of Hamil¬ 
tonian (2) [41]. For our exponentially-decaying interaction 
potential (9), we have £[wexp](—7") = — T)~^. In 

the LL limit, £[z(;](—T) = |1, and the term proportional to 
£[u'exp](—in Eq. (A.ll) is zero. Taking the derivative of 
Eq. (A. 7), now with respect to Y\ yields the equation for the 
optimal Y*, 

27r5(0)y* = ^(^liTlT-) (A.12) 

and hence W* and V*. Erom these, Q and R are evolved 
forward in time by a step dr. 

As mentioned above, we fix the gauge of IT*) such that I = 
1. This is achieved by choosing an arbitrary R and an anti- 
hermitian K, and setting Q = K — ^R^R. It is then possible 
to construct an update for Q and R that preserves this gauge 
exactly, 

R{t + dr) = R{t) — dr W* (r) 

K{t + dr) = K{t) Y y (i?(T)^fF*(r) - W*{t)'^R{t)'^ . 

(A.13) 

Indeed one can check that Q{T + dT) = Kir + dr) — 


dr)^R {t + dr). To first order in dr equals the update R{t + 
dr) = R{t) — dTW*{T) and Q{t + dr) = Q{t) — dTV*{T) 
(remember Z = 1 ). 

Eor the homogeneous case, it is also possible to manually 
regauge the matrices Q, R after each update step. Suppose 
Q, R are ungauged CMPS matrices. Regauging them such 
that Z = 1 is done by first calculating the left eigenvector 
(Z| of T to the eigenvalue a with largest real part. The state 
jT*) is then renormalized hy Q ^ Q = Q — ^1. R and 
Q are then transformed according to Qi — and 

Ri = Z^/^i?Z“^/^. It is easy to check that for these matrices, 
the left eigenvector Z with eigenvalue 0 is indeed the identity 
matrix. 

The inverse (—T)p^ in Eq. (A.ll) is to be understood as a 
pseudo-inverse acting on the vector space orthogonal to |r)(Z|. 
This regularizes the infinite energy content coming from sum¬ 
ming up all local energy contributions to the left of a partic¬ 
ular position X. The action (Zc|(—T)p^ = {k\ on an arbitrary 
given vector {k\ is computed from solving the inhomogeneous 
linear equation system {k\ [1 — |r)(Z|] = —{k\ [T — |r)(Z|]for 
{k\, using a sparse solver [42] for non-hermitian equation sys¬ 
tem. A similar approach without any pseudo-inverse is used 
to calculate {rj — T)~^. To get a stable update we have used an 
implicit Euler scheme to update Q and R. The overall oper¬ 
ational cost of the procedure explained above is 0{D^). De¬ 
pending on the parameters rj and g, time steps dr have to be 
chosen as small as dr = 10 “^. 
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